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ABSTRACT 

Petrosian magnitudes were designed to help with the difficult task of determining a galaxy's total 
light. Although these magnitudes (taken here as the flux within 2Rp, with the inverted Petrosian index 
l/r)(Rp) — 0.2) can represent most of an object's flux, they do of course miss the light outside of the 
Petrosian aperture (2Rp). The size of this flux deficit varies monotonically with the shape of a galaxy's 
light-profile, i.e., its concentration. In the case of a de Vaucouleurs R 1 ' 4 profile, the deficit is 0.20 
mag; for an i? 1 / 8 profile this figure rises to 0.50 mag. Here we provide a simple method for recovering 
total (Sersic) magnitudes from Petrosian magnitudes using only the galaxy concentration (Rqq/R^q or 
i?8o/i?2o) within the Petrosian aperture. The corrections hold to the extent that Sersic's model provides 
a good description of a galaxy's luminosity profile. We show how the concentration can also be used 
to convert Petrosian radii into effective half-light radii, enabling a robust measure of the mean effective 
surface brightness. Our technique is applied to the SDSS DR2 Petrosian parameters, yielding good 
agreement with the total magnitudes, effective radii, and mean effective surface brightnesses obtained 
from the NYU-VAGC Sersic R}l n fits by Blanton et al. (2005). Although the corrective procedure 
described here is specifically applicable to the SDSS DR2 and DR3, it is generally applicable to all 
imaging data where any Petrosian index and concentration can be constructed. 

Subject headings: galaxies: fundamental parameters — galaxies: structure — methods: analytical — 
methods: data analysis 



1. INTRODUCTION 

Galaxies are known to possess well-defined correlations 
between their various structural and kinematic parameters 
(e.g., Faber & Jackson 1976; Tully & Fisher 1977; Djorgov- 
ski & Davis 1987; Dressier et al. 1987; Caon, Capaccioli, 
& D'Onofrio 1993; Graham, Trujillo, & Caon 2001; De Ri- 
jcke et al. 2005; Matkovic & Guzman 2005). These empir- 
ical "scaling-laws" provide key observational constraints 
needed to test current theoretical models of galaxy for- 
mation and evolution. Obviously even the most basic of 
these, the luminosity-size relation (e.g., Dutton et al. 2005; 
Mcintosh et al. 2005), relies on our ability to accurately 
measure robust photometric parameters. 

In this vein, the need for a more unified approach to 
galaxy photometry has recently been highlighted by Cross 
et al. (2004). They noted that much of the discrepancy 
in galaxy magnitudes (and sizes) between various groups 
is because of the varying methodology applied. For exam- 
ple, some Authors use Kron magnitudes, others Petrosian 
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magnitudes, some extrapolate fitted models to large radii, 
while others use somewhat limited aperture photometry. 
This can impact significantly on global measures of the 
galaxy population, such as the luminosity function (e.g., 
Norberg et al. 2002; Blanton et al. 2003a; Driver et al. 
2005), the color-magnitude relation (e.g., Scodeggio 2001; 
Chang 2005, and references therein), and the luminosity 
density (e.g., Yasuda et al. 2001; Cross et al. 2001). It also, 
for example, impacts on studies of the supermassive black 
hole mass function derived using the galaxy luminosity- 
black hole mass relation (e.g., McLure & Dunlop 2004, 
Shankar et al. 2004, and references therein). Perhaps less 
obvious however, are the consequences for the calcula- 
tion of size and surface brightness distributions (e.g., Kor- 
mendy 1977; Cross et al. 2001; Shen et al. 2003; Driver et 
al. 2005). These are typically derived from the half-light 
radius which in turn depends critically on an accurate as- 
sessment of the total flux. If the magnitude is underes- 
timated, the size and surface brightness distributions will 
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be affected. 

One of the great strengths of the Petrosian (1976) index, 
the average intensity within some projected radius divided 
by the intensity at that radius, and similarly the radii 
themselves corresponding to some fixed Petrosian index, 
is that they do not depend on a galaxy's distance. That 
is, because surface brightness dimming does not change 
the shape of a galaxy's light-profile, it does not affect the 
Petrosian index, nor does it affect the observed galaxy 
concentration. Furthermore, due to the Petrosian index's 
ability to define aperture sizes which contain the bulk of 
an object's light, and due to the large influx of small, faint 
images of high-redshift galaxies that are now available, 
the Petrosian index has experienced a resurgence (e.g., 
Wirth, Koo, & Kron 1994; Bershady, Lowenthal, & Koo 
1998; Dalcanton 1998; Takamiya 1999; Bershady et al. 
2000; Volontcri, Saracco, & Chincarini 2000; Blanton et 
al. 2001; Lubin & Sandage 2001; Consclicc, Gallagher, & 
Wyse 2002; Yagi et al. 2002). Strauss et al. (2002; their 
section 3.2) do however stress the fact that a different 
fraction of galaxy light is missed depending on whether 
a galaxy has an i? 1 / 4 light-profile or an exponential light- 
profile, and they emphasize the subsequent need to ac- 
count for this in analyses of the Sloan Digital Sky Survey 
{SDSS; York et al. 2000) galaxy data. 

In an effort to account for the flux missed by Petrosian 
apertures, this paper outlines a corrective procedure to 
convert Petrosian magnitudes into total (Sersic) galaxy 
magnitudes, as provided by popular codes such as GIM2D 
(Simard 2002), GALFIT (Peng et al. 2002), and BUDDA 
(de Souza, Gadotti, & dos Anjos 2004). The key to doing 
this lies in the 'shape' of a galaxy's stellar distribution, 
that is, its concentration. This quantity may be obtained 
with or without the use of a fitted light-profile model. 

In the absence of measurement errors, all i? 1 / 4 
light-profiles have exactly the same concentration index 
-R90 /R50 — the ratio of radii containing 90% and 50% of 
the Petrosian flux (e.g., Blanton et al. 2001; Strauss et 
al. 2002; Goto et al. 2003). The same is true for an ex- 
ponential i? 1 / 1 light-profile, although the specific value of 
the concentration will be different in this case. It follows 
that the observed range of galaxy concentrations, if not 
due to errors, reflects a range of light -profile shapes; that 
is, galaxies do not simply have exponential or i? 1//4 light- 
profiles. 

Recognizing this in the SDSS data, Blanton et al. 
(2003b) adopted Sersic's (1963, 1968) i? 1 /™ model 3 to rep- 
resent the range of galaxy light-profile shapes and pro- 
vide estimates of their total luminosities, sizes, and surface 
brightnesses. Indeed, in the case of (dwarf and ordinary) 
elliptical galaxies, such an approach is crucial if one is 
to properly understand the various relationships between 
such terms (e.g., Graham & Guzman 2003, their Section 
4). However, for low signal-to-noise data or where the 
spatial resolution is lacking, it can become difficult to ob- 
tain reliable Sersic fits. One therefore needs an alterna- 
tive strategy to obtain the total galaxy flux and associated 
half-light terms. 

Using only the observed concentration, Rqq/R^q, within 
the Petrosian aperture, we provide an easy prescription 



to recover the total (Sersic) flux from the Petrosian flux 
while maintaining the distance-independent qualities of 
the Petrosian system. We also explain how one can re- 
cover the effective radii and associated surface brightness 
terms. The corrective formula presented here are not only 
valid for pure i? 1 / 4 or exponential profiles, but applicable 
to galaxies having intermediary light-profile shapes and a 
range of more extreme stellar distributions. 

2. PETROSIAN RADII AND MAGNITUDES 

The Petrosian (1976) index was initially introduced with 
the goal of measuring galaxy evolution. It gained addi- 
tional popularity from its potential to determine the cos- 
mological parameters (e.g., Djorgovski & Spinrad 1981; 
Sandage & Perelmuter 1990). Indeed, under the assump- 
tion of structural homology, it provided a means to obtain 
"standard rods" that could be used to constrain cosmogo- 
nic models. Nowadays it is often used as a tool for defining 
aperture sizes from which to measure galaxy magnitudes. 

The Petrosian index, rj{R), is a function of a galaxy's 
projected radius R, and can be written as 



v(R) 



2 J*I(R')R'dR' _ L(<R) _ (I)! 



(1) 



R 2 I(R) nR 2 I(R) I(R) ' 

where I(R) is an object's (projected) intensity at some 
radius R, and (I)r is the average intensity within that ra- 
dius. Following, for example, Blanton et al. (2003b), who 
modeled 183,487 SDSS galaxies, we adopt Sersic's (1963) 
R}l n model to represent the possible range of light-profiles 
I(R). This can be written as 

V 

-b. 



I(R) = I c exp 



^ -1 



(2) 



where I e is the intensity of the light profile at the effec- 
tive radius R c , and n defines the 'shape' of the profile. The 
term b n is simply a function of n and chosen 4 to ensure 
the radius R e encloses half of the profile's total luminos- 
ity. Using the substitution x = b n (R/R c ) 1 ^ n , and thus 
R = x n R e /(b n ) n and dR = {R e n/{b n ) n )x n - 1 dx, the Pet- 
rosian index reduces to 

2nj(2n, x) 



r](x,n) 



(3) 



where 7(2n,x) is the incomplete gamma function defined 

as 

"Mt. (4) 



■y(2n,x) = / 
Jo 
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The inverted Petrosian index, l/r](R), is more com- 
monly used in the literature (e.g., Bershady, Jangren, & 
Consclicc 2000; Blanton et al. 2001) and is shown in Fig- 
ure 1. It has a value of 1 at R = and falls to zero at 
large radii. Although the Petrosian index can be used to 
determine a galaxy's 'Petrosian radius' Rp — the radius 
where the index equals some fixed value — it is important 
to determine the relation between this radius and the ef- 
fective radius R c : as shown inset to Fig. 1, Rp/Re varies 
with Sersic index n for a fixed 1/rj. Constant multiples of 
Rp have been used as a means to define an appropriate 



aperture size for purposes of deriving galaxy magnitudes. 

3 A useful compilation of various Sersic expressions can be found in Graham & Driver (2005). 

4 The value b n is such that T(2n) = 27(2n, b„), with 7 and T the incomplete and complete gamma functions respectively (Ciotti 1991). 
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For example, some Authors have chosen to use 2Rp with 
l/r)(Rp) = 0.2 (e.g., Bershady et al. 2000; Blanton et al. 
2001) and others 3Rp with l/ry(i?p) = 0.5 (e.g., Conselice 
et al. 2002, 2003). The SDSS consortium adopted the for- 
mer criteria. From Figure 1, one can see that different 
light-profile shapes reach constant values of 1/n at differ- 
ent fractions of their effective half light radii R c . (The 
absolute value of R e , and I c , are not important, only the 
ratio R/R c and the value of n determine 1/rj, see equa- 
tion 3.) 

This is shown more clearly in the inset Figure where 
one can see, as a function of profile shape n, the number 
of effective radii that 1/n = 0.2 and 0.5 correspond to. If 
n=4, then l/f](Rp) = 0.2 occurs at R P = 1.82i? c . How- 
ever, when using l/rj(Rp) = 0.5, if n=A then Rp equals 
only 0.16i? e 5 . Although this particular radius is multi- 
plied by a factor of three before determining the Petrosian 
magnitude, this still amounts to an aperture of size less 
than 0.5i? c - This particular Petrosian magnitude there- 
fore greatly underestimates an object's total magnitude. 

For a Sersic profile, the Petrosian magnitude, mp, is 
given by the expression 

r e fc„ 

m(< NRp) = [i c — 51ogi? e — 2. 5 log 2-7m 

L \yn) 



-7(2n,x P ) 



results in magnitude differences of 0.5 mag when n ~ 2.5; 
1.25 mag when n = 4; 2 mag when n ~ 5.5 and con- 
siderably worse for galaxies with yet higher values of n. 
Obviously such an approach to determine galaxy magni- 
tudes should be used with caution. When dealing with 
dwarf galaxies, because of their faint central surface bright- 
nesses, the sky flux can often dominate at the radius 
where l/r)(R P ) = 0.2. The use of l/r}(R P ) = 0.5 is thus 
more practical, and for galaxies with n < 2 the bulk of 
their flux is still recovered. Overall, however, the use of 
l/rj{Rp) = 0.2 does much better at recovering a galaxy's 
true magnitude 6 . When n=4, the magnitude difference 7 
is only 0.20 mag, rising to 0.64 mag when n = 10. Things 
arc better for profiles with smaller values of n and if one 
considers the profiles truncate at 5i? e (Figure 2b). 

3. RECOVERING EFFECTIVE RADII AND TOTAL (SERSIC) 
MAGNITUDES 

In this section we will only consider the case where 
l/n(Rp) = 0.2, using 2i?p to define the Petrosian mag- 
nitude. The merits of this precise definition, which was 
used by the SDSS team, are exalted in Strauss et al. (2002, 
their section 3.2). 



(5) 

where xp = 6 n (A r i? P /i? e ) 1 /™, N is the multiplicative fac- 
tor (usually 2 or 3), and /x e = — 2.5 log 7 . The total mag- 
nitude is obtained by replacing 7(2n, xp) with the (com- 
plete) gamma function T(2n). The extent to which the 
Petrosian magnitude underestimates the total magnitude 
is shown in Figure 2 under two conditions. The first is 
when the Sersic profiles extend to infinity (panel a). Al- 
though galaxies are recognized not to have sharp edges, the 
second condition assumes that the Sersic profiles truncate 
at 5i? e (panel b). There is no physical justification for a 
truncation at 5i? c ; this is simply chosen in order to explore 
the magnitude deficit if the light-profiles don't extend to 
infinity. While elliptical galaxy light-profiles continue into 
the background sky-noise (e.g., Caon et al. 1993), there is 
evidence that some disk galaxies may truncate at ^4 scale- 
lengths, or change their exponential slope at these radii 
(van der Kruit 2001; Pohlcn et al. 2004; Erwin, Beckman, 

6 Pohlen 2005), but see Narayan & Jog (2003) and also 
Bland-Hawthorn et al. (2005) who present a light-profile 
for the late-type disk galaxy NGC 300 which extends to 
10 scale-lengths. Figure 2 thus provides a boundary of 
sorts to the extent that Petrosian magnitudes may under- 
estimate a galaxy's total magnitude as a function of the 
underlying profile shape. To avoid possible confusion, we 
note that R e shall always refer to the value associated with 
the R}l n model extended to infinity. Thus, when the pro- 
file is assumed to truncate at 5i? e , the value of R c does 
not change. 

One can see from Figure 2a that the use of l/77(i?p)=0.5 

5 Although Section 3.3 of Bershady et al. (2000) reports that l/r](Rp) = 0.5 roughly corresponds to Rp = IRe, this is only true for Sersic 
profiles with n < 2 (see Figure 1). Exact values for Rp, in terms of R e , when l/rj(Rp) = 0.5 are given in Bershady et al. (1998) for n = 0.5, 1, 
and 4. 

° The use of l/r?(.Rp) = 0.2 has also been adopted by some Authors because it results in a minimal variation of Rp/R c with n. For n < 4, this 
value is around 2, and so to approximate the total light, one could measure the light within Rp/2 and multiply by 2. 

7 Bershady et al. (2000) reported a difference of 0.13 mag, but this is appropriate for an n = 3 profile rather than an n = 4 profile. The value 
of 0.13 mag was based on photometry of IRAF artdata simulations, and therefore is likely to be a round-off or discretization error in the IRAF 
rendering. 

8 If an n = 1 and n = 4 profile arc integrated to infinity, rather than only 2_Rp, then Rgo = 2.32i? c and 5.55-R e respectively. 



3.1. Concentration 

If one had some way of knowing the underlying light- 
profile shape 'n', then from Figure 2 one could correct the 
Petrosian magnitudes for the missing flux beyond 2Rp. 
Given that the Petrosian index was developed in part to 
avoid fitting a model to an object's light-profile, one may 
prefer not to fit the R}l n model, or indeed one may not 
have the resolution to do so. Conveniently, central light 
'concentration' is monotonically related to the shape of a 
light-profile (Trujillo, Graham, & Caon 2001a; Graham et 
al. 2001). One can therefore use concentration as a proxy 
for the value of n. 

The SDSS consortium have been using a ratio of two 
radii (i?g and i? 50 ) as a measure of an object's concen- 
tration, both of which are available from the SDSS public 
data releases. These radii enclose 90% and 50% of the Pet- 
rosian flux and their ratio is shown in Figure 3 and Table 1 
as a function of the underlying light profile shape n. The 
use of such radii avoids the inner seeing affected part of 
a light-profile and also avoids the outer noisier part of a 
profile while still providing a useful range of concentrations 
for the various galaxies. For an n = 1 and n = 4 profile 8 , 
R90/R50 equals 2.29 and 3.42 respectively. The inverse 
ratio is sometimes used, giving 0.44 and 0.29 (Blanton et 
al. 2001, their section 4.5). Not surprisingly, the R 9n /R 50 
concentration index correlates with galaxy type and also 
color (e.g., Strateva et al. 2001; Kauffmann et al. 2003). 

3.2. Magnitudes 
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Combining Figures 2 and 3, Figure 4a shows the re- 
quired magnitude correction in order to account for the 
missing flux beyond the Petrosian aperture 2Pp, when 
l/j](Rp) = 0.2, as a function of concentration within the 
Petrosian aperture (i.e., curve (hi) from Figure 3). When 
n = 4, i? 90 /i?5o = 3.42, and Amag = 0.20 and 0.07 mag 
to recover the total flux and that within 5P C , respectively 
(see Table 1). 

To an accuracy of ~ 0.01 mag, over the Sersic interval 
0.1 < n < 10, the missing flux (Figure 4a) can be approx- 
imated by the expression 



Am w Pi exp (P 90 /P 5 o) 



N-P2 



to p - m x , 



(6) 



where Pi and P2 equal 5.1 xl0~ 4 and 1.451 in order to re- 
cover the total flux (Figure 5a), and 5.9x 10~ 5 and 1.597 in 
order to recover the flux within 5P C (Figure 5b). Here, top 
is the Petrosian magnitude and m x is the corrected mag- 
nitude which, to an accuracy of ~0.01 mag, is equivalent 
to the Sersic magnitude (tos). 

3.3. Radii 

The effective radii containing half of the total (Sersic) 
flux can be computed in a number of ways. 

Figure 4b is the combination of Figure 3 and Figure l's 
inset figure, and allows one to determine the effective 
radius R c from the Petrosian radius and concentration 
P90/P50 within the Petrosian aperture. Alternatively, 
from the total galaxy magnitude (approximated by m x ) 
one may empirically determine the aperture containing 
half a galaxy's light and therefore obtain the half-light 
radius this way. 

Alternatively still, for 0.1 < n < 10, the Pso/Pe curve 
seen in Figure 4c can be approximated in terms of Pet- 
rosian concentration by the expression 

R5 ° ~R X , (7) 



Pc 



P3 (P90/P50) 



Pi 



where P3 and P4 equal 8.0 x 10 6 and 8.47, respectively. 
This approximation is shown in Figure 5c. 

3.4. Surface brightnesses 

It turns out one can also easily transform ^ 50 , the sur- 
face brightness at P50, and (^)5o, the mean surface bright- 
ness within P50, into the effective surface brightness, /i , 
and the mean effective surface brightness (n) c - From 
Graham & Driver (2005, their equation 6 and their sec- 
tion 2.2), one has that 

l/n 



M50 — Mo 



2.5b» 
ln(10) 



and 



(M)50 - (m)o = 2.5 log 



P30 

Re 
Pc 



1 



7(2n,x c ) 



(8) 



(9) 



7(2n,x 5 o) 

where x$o = b n {Rbo/ Rc) 1 ^ 11 and x c = b n . One can see that 
for a given value of n, the ratio Pso/P c (Figure 4c) is suf- 
ficient to solve equations 8 and 9. Again using the relation 
between n and P90/P50 (Figure 3, curve (hi)), Figures 4d 

9 SDSS-DR2: http://www.sdss.org/dr2/ 

10 NYU-VAGC: http://sdss.physics.nyu.edu/vagc/ 



and 4e show the surface brightness differences as a func- 
tion of the concentration Pgo/Pso- In practice, one can 
directly measure the surface brightness at P e and/or the 
mean surface brightness within P , the latter of which can 
also be computed from the expression L to t = 27rp2( /) e . 
Approximations to equations 8 and 9 are therefore not 
given. 

Because the above values of P c ,/i c , and (^} c are those 
pertaining to a non-truncated Sersic profile, equations 7-9 
are independent of any possible profile truncation beyond 
the Petrosian aperture. 

4. APPLICATION 

4.1. Demonstration with SDSS data 

The SDSS consortium adopted a slightly modified form 
of the Petrosian index. Rather than dividing the average 
intensity within Pp by the intensity at Pp (Equation 1), 
they divided the average intensity within Pp by the aver- 
age intensity within 0.8-1.25Pp. This was done so as to 
reduce the sensitivity of the index to noise and (possible) 
real small-scale fluctuations in the light -profile (Strauss 
et al. 2002). As a result, their definition of the Petrosian 
index is given by the expression 



rj(x,n) 



7(2n, XL25R) - j(2n, x Q . 8R ) 
[(1.25) 2 - (0.8) 2 ] 7 (2n,a;) ' 



(10) 



where xi. 2 5_r = 6„(1.25P/P c ) 1/n and x . 8 r = 
6„(0.8P/P C ) 1 /™. 

In general, where the SDSS Petrosian index (Equa- 
tion 10) equals 1/0.2, the associated Petrosian radii are 
slightly smaller than the radii obtained previously with 
Equation 3 equal to 1/0.2 (see Table 1). The required 
magnitude corrections are therefore slightly greater when 
using the SDSS definition. The overall appearance of Fig- 
ures 1-4 do not change, but the exact numbers are different 
by a few percent and therefore provided in the middle sec- 
tion of Table 1 . The approximations given by equations 6 
and 7 also require a slight modification. Using the SDSS 
definition of the Petrosian index (equation 10), to recover 
the total magnitude one now has Pi = 4.2 x 10~ 4 and 
P2 = 1.514 (Figure 5d), and to recover the flux within 
5P e Pi = 8.0 x 10~ 5 and P 2 = 1.619 (Figure 5e). The 
values for P3 and P4 are 6.0 x 10~ 6 and 8.92, respectively 
(see Figure 5f). 

We have tested the applicability and need for these cor- 
rections using real data from the SDSS consortium. Specif- 
ically, we have taken the Petrosian magnitudes, 50-percent 
Petrosian radii, and surface brightnesses (top, P50, (m)so) 
from the SDSS Second Data Release 9 (Abazajian et al. 
2004) and the Sersic-derived total magnitudes, effective 
radii, and surface brightnesses (tos, R e , (A*)e) from the New 
York University Value-Added Galaxy Catalog (NYU- 
VAGC, Blanton et al. 2005) 10 . From this sample we then 
removed those galaxies whose Petrosian derived P50 val- 
ues are biased high by seeing; we crudely did this by ex- 
cluding galaxies with P c < 5" (90% of SDSS DR2 see- 
ing was between 1.0-1.6"). In deriving Sersic parameters, 
Blanton et al. (2005) prevented the fits from obtaining 
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Sersic indices greater than 5.9 (and less than 0.2). To 
avoid any objects piled up at this boundary and which 
therefore may not have accurate Sersic quantities, we only 
use galaxies with n < 5.8 (and greater than 0.21). This 
left us with a sample of 16128 galaxies, more than enough 
to test our method 11 . 

Figure 6 plots the difference between the Petrosian and 
Sersic values as a function of the (SDSS-tabulated) con- 
centration indices -R90/-R50 within the Petrosian apertures. 
A similar figure using Sersic n rather than concentration is 
given in Blanton et al. (2003a, their figure 14). The solid 
curves show the expected differences that we have derived. 
The dashed curves show the selection boundary imposed 
upon the data due to the restriction used by Blanton et 
al. (2005). Confining n to values smaller than 5.8 can pre- 
vent an accurate recovery of the total Sersic flux in large 
galaxies, and (artificially) prevents Sersic magnitudes get- 
ting more than 0.37 mag brighter than the SDSS-defincd 
Petrosian magnitudes (see Table 1, middle section). Of 
course, measurement errors in either the Pertrosian or 
Sersic magnitude, and also real deviations from a Sersic 
profile, can result in differences outside of this selection 
boundary. 

In Figure 6a, the Sersic magnitudes are clearly brighter 
than the Petrosian magnitudes for concentrations greater 
than about 2.7 (n ~ 1.7). In Figure 6b it is obvious 
that our corrected Petrosian magnitudes agree with the 
Sersic magnitudes, at least to a concentration of ~3.3 
(n <~ 4) at which point the selection boundary from the 
Sersic fits start to dominate the figure. In the concentra- 
tion bin 2.75 < i?9o/i?5o < 3.0, the mean difference be- 
tween the Petrosian and Sersic magnitudes, (mp — ms), 
equals 0.087±0.005 mag: close to the expected offset 
of 0.06 mag at -R90/-R50 = 2.87 (the mean concentra- 
tion inside this bin). In this same bin, (m x — ms) = 
0.027 ±0.005 mag. In the next bin, 3.0 < R90/R50 < 3.25, 
(mp— ms) = 0.144±0.004 mag, with the expected value at 
-R90/-R50 = 3.12 (the mean concentration inside this bin) 
equal to 0.13 mag. Similarly, (m x - m s ) = 0.028 ± 0.004 
mag in this bin. The mean value of i?go/i?5o in the bin 
3.25 < i?9o/i?5o < 3.5 is 3.34, for which the expected value 
of (mp — ms) is 0.22 mag. However, the measured value 
is only 0.148±0.006 mag. The reason for this is partly 
due to the selection boundary but primarily due to the 
underestimation of the Sersic magnitude by Blanton et al. 
(2005). Figure 9 from Blanton et al. reveals that at n = 4, 
corresponding to R 9n /R 50 — 3.35, their Sersic fluxes are 
underestimated by ~7%. This therefore accounts for the 
value of 0.148±0.006 instead of ~0.22, and for the negative 
value of (m x — ms) = —0.063 ± 0.006 mag in this bin. 

For a system with n = 5, Blanton et al. (2005) typically 
obtained n =4.0-4.6 and recovered only ~90% of the ac- 
tual flux from their simulated test galaxies. The situation 
is systematically worse for galaxies with higher values of 
n and accounts for much of the apparent over-correction 
in Figure 6b at high-concentrations. Due to this problem, 
the simple corrective procedure presented here to allow for 
the missing flux outside of the Petrosian aperture actu- 
ally provides a means to acquire (total) Sersic magnitudes 

11 For those who need or wish to use the full galaxy sample, the seeing-i 
and 2001c for a prescription to do this). The Sersic radii R e were obta 
Sersic profiles and therefore need no further seeing-correction. 



which are more accurate than those obtained from the di- 
rect application of Sersic models to the data. For a typical 
elliptical galaxy, our method improves the accuracy on the 
galaxy magnitude from a couple of tenths of a mag to a 
couple of hundredths of a mag (assuming that the Sersic 
profile within the Petrosian aperture continues outside of 
it). 

We do note that >99% of the SDSS Main Galaxy Sample 
(m r < 17.77) have R90/R50 < 3.5 (Blanton et al. 2003b; 
Nakamura et al. 2003). Thus, cases where corrections are 
most severe are rare. This, however, is not to say that the 
high-concentration objects are uninteresting. Due to the 
Mbh-concentration relation (Graham et al. 2001), they are 
the galaxies expected to have the most massive supermas- 
sive black holes. 

In passing we also note that Figures 6a and 6b reveal 
that the SDSS Petrosian magnitudes are too faint, in the 
sense that they underestimate the total galaxy flux, at 
the high- luminosity end. Correcting for this increases the 
galaxy number counts at the bright-end of the (Petrosian- 
derived) luminosity function (see Blanton et al. 2003a). 

Figure 6c shows the difference between the Petrosian 
radii R$o and the Sersic effective radii. There is a clear 
offset at the high-concentration end of the diagram, al- 
beit constrained within the selection boundary which pre- 
vented Sersic indices greater than 5.8 and thus, from Ta- 
ble 1, i?5o/i? c < 0.53. At the concentration R 90 /R 5n = 
3.4, n ~ 4.5 and from Figure 9 in Blanton et al. (2005) 
one can see that their Sersic effective radii are only ^80% 
of the true effective radii. Allowing for this, the average 
ratio R x /R c should be ~1.25 at -R90/-R50 = 3.4, exactly 
as observed in Figure 6d. That is to say, the apparent 
mis-match at high-concentrations is understood: it arises 
from the systematics which Blanton et al. (2005) identified 
in their fitted Sersic quantities. 

Figure 6e reveals that the early-type galaxies, defined 
to be those with concentrations greater than 2.86 (e.g., 
Shimasaku 2001; Nakamura et al. 2003), have Petrosian 
surface brightnesses which are clearly brighter than the 
Sersic-derived surface brightness. Figure 6f shows that 
the procedure does a good job at correcting for the miss- 
ing flux beyond the Petrosian apertures, at least until a 
concentration of around 3.3 where the selection boundary 
and the underestimation of the Sersic flux once more come 
into play. 

4.2. Application in CAS space 

Another popular setup is that used in the CAS system 
(Conselice 2003), which measures the concentration (C), 
asymmetry (A), and stellar/star-forming dumpiness (S) of 
a galaxy's light distribution. Within the CAS code, Pet- 
rosian apertures have sizes of 1.5i?p, with l/ry(i?p) = 0.2. 
Concentration is also defined slightly differently from that 
used by the SDSS consortium: it is the ratio of radii con- 
taining 80% and 20% (rather than 90% and 50%) of the 
flux within the Petrosian aperture (Bershady et al. 2000). 

Theoretically, it doesn't matter what choice of radii one 
uses to define concentration. For example, radii contain- 
ing 100% and 50% of the flux within the Petrosian aper- 
ture would work. In practice, the use of a large upper- 

icted radii R50 will first need to be corrected (see Trujillo et al. 2001b 
d by Blanton et al. (2005) from the application of seeing-convolvcd 
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percentage and a small lower-percentage leads to a greater 
range of concentrations and thus a clearer distinction be- 
tween different profile types. However, not using too large 
an upper-percentage allows one to work with a less noisy 
part of the light-profile. Not using too small a lower- 
percentage means that one is less vulnerable to the effects 
of seeing. An analysis of the optimal concentration index 
(e.g., Graham et al. 2001) is, however, beyond the intended 
scope of this paper and will be addressed elsewhere. We 
note that such a search need not be limited to differing per- 
centages, but could include radii based on different values 
of the Petrosian index. Here we simply provide the cor- 
rections using the popular SDSS and CAS definitions for 
the concentration index and Petrosian index/ aperture. 

The various size and flux transformation and correc- 
tive terms for the specific Petrosian setup (1.5i?p, with 
l/rj(Rp) = 0.2) and concentration (-R80/-R20) used by the 
CAS code of Conselice et al. (2003) are given in the lower 
section of Table 1. To approximate the total Sersic mag- 
nitude, over a range in n from 0.1 to 10, this missing flux 
Am can be approximated by 

Am w -0.096 + 0.021(i? 80 /i? 20 ) + 0.00 44(i? 80 /i? 20 ) 2 
= m P -m x . (11) 

To obtain the flux within 5R e , the following provides a 
good approximation 

Am « -0.012 - 0.0013(i? 80 /i?2o) + 0.00346( J R 80 /i? 2 o) 2 - 

(12) 

The effective Sersic radius R c is given by 
R 50 /R e « 0.903 + 0.102(R 80 /R 20 ) - 0.0276(i? 80 /i? 20 ) 2 

+0.00119(i? 80 /i? 20 )fl3) 

4.3. General applicability 

Caon et al. (1993) have shown Sersic's R}l n model fits 
early-type galaxies remarkably well (see their figure 2) 
down to faint surface brightness levels (p,B ~ 27 mag 
arcsec -2 ). Spiral galaxies, however, usually have two 
clearly distinct components, namely a bulge and a disk. It 
is therefore pertinent to inquire how the above corrections 
may apply when a galaxy is clearly better represented by 
two components, as is the case for the intermediate-type 
galaxies. 

When dealing with luminous galaxies, the rough divide 
between early-type and late-type galaxies has been taken 
to occur at R90/R50 = 2.86 (n ~ 2) (Shimasaku 2001; 
Nakamura et al. 2003; Shcn et al. 2003, using the SDSS 
definition of Petrosian index and aperture). 12 It should 
be kept in mind that this concentration based definition of 
galaxy type is only applicable to luminous systems. Dwarf 
elliptical galaxies have exponential-like profiles and thus 
low concentrations (e.g., Graham et al. 2001; their figure 
8), but are obviously not late-type galaxies. The magni- 
tude correction for galaxies having concentrations of this 
size or smaller is less than 0.06 mag (Table 1). In late-type 
galaxies, the exponential disk, rather than the bulge, dom- 
inates the flux in the outer parts. Whether or not these 
disks continue for many scale-lengths or truncate at a few 
scale-lengths, if one applies the corrections presented here, 
one will not over-correct the galaxy flux by more than 0.06 

12 Other Authors have used n = 2.5 (e.g., Blanton ct al. 2003b) or 
the rough divide between luminous early- and late-type galaxies. 



mag for the 2-componcnt late-type galaxies, nor overesti- 
mate the half-light radius by more than 7% (Table 1). If, 
however, one is able to identify the disk-dominated galax- 
ies from their blue rcstframe colors, or indeed from their 
low concentrations (assuming no low luminosity elliptical 
galaxies reside in one's sample, or folding in asymmetry 
to separate the galaxy classes) , then one may instead elect 
to apply the small corrections applicable for systems with 
an n=l (exponential disk) light-profile. Given the cor- 
rections for such stellar distributions are small, a simple 
answer may then be to only apply the various corrections 
to the red galaxy population, or to systems with (SDSS) 
concentrations greater than say 2.86. 

We caution that intermediate objects, such as lenticu- 
lar galaxies, with half their light coming from their disk 
and half from their bulge may not be well approximated 
with a single Sersic profile. In such cases the concentra- 
tion prescription given here should be used with care as it 
is intended for systems which can be approximated with 
a single Sersic function. In general, however, Section 4.1 
and in particular Figure 6 reveal that the method devel- 
oped here seems to work rather well. 

A different issue pertains to how one actually mea- 
sures the Petrosian index. For a disk system, the Pet- 
rosian index is independent of disk inclination if one 
uses appropriate elliptical apertures reflecting the incli- 
nation of the disk. If one instead measures the index 
from the slope of the major-axis light-profile, such that 
2/77 = d[lnL(R)]/d[\nR] (Gunn & Oke 1975), then again 
the index is not dependent on the inclination of the disk. 
However, the use of circular apertures applied to highly- 
inclined disk galaxies will result in erroneously high con- 
centrations compared to what these values would be if 
such galaxies were viewed with a face-on orientation. Al- 
though we have not quantified this effect, comparisons be- 
tween eye-ball morphology and concentration index yields 
reasonable agreement for the late-type galaxies (e.g., Shi- 
masaku et al. 2001), suggesting that in practice this effect 
is not a significant problem. 

Lastly, we note that if image distortion due to seeing is 
not corrected for and is such that the PSF has altered the 
Petrosian 50% radius R50 and thus -R90/-R50 (or i? 8 o/i?2o) 
to the extent that the measured value no longer reflects the 
intrinsic galaxy concentration, then obviously one should 
not be including such objects in plots or analyses that 
use this quantity. The concentration-based corrections 
outlined in this study cannot be applied in such circum- 
stances, at least until the radii are corrected for seeing 
(e.g., Trujillo ct al. 2001b, c) We also direct readers to 
Blanton et al. (2001) who present a study of how the SDSS 
Petrosian magnitudes are affected by seeing, and Wirth et 
al. (1994) who considered the affects of both poor image 
resolution and low signal-to-noise on image concentration. 

5. SUMMARY 

In commenting on the adopted Petrosian system of the 
SDSS consortium, Strauss et al. (2002) wrote: "Any sci- 
entific analysis that uses the [SDSS] redshift survey must 
consider how the fact that a different fraction of light is in- 
cluded for elliptical and spiral galaxies affects the result" . 

"?50 = 2.6 (e.g., Strateva et al. 2001; Kauffman ct al. 2003) to mark 
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This was both a recognition that a) their particular Pet- 
rosian system likely captures a different fraction of galaxy 
light compared to other surveys and that b) any practi- 
cal definition of Petrosian parameters will inevitably miss 
some fraction of a galaxy's light, and a different fraction 
in different galaxies. For example, a galaxy with an expo- 
nential light-profile will have 99% of its total flux encom- 
passed by the SDSS Petrosian magnitude, but if a galaxy 
has an i? 1 / 4 profile, then the flux beyond the SDSS Pet- 
rosian aperture accounts for an additional 0.22 mag, i.e., 
roughly 1/5 th of that galaxy's total light. In this paper 
we outline a simple method for accurately recovering the 
missed flux. 

Differences between galaxy light-profile shapes result in 
a range of galaxy concentrations. We have shown how 
these can be used to correct for the flux outside of the Pet- 
rosian apertures and thus provide a galaxy's total (Sersic) 
magnitude. The corrections presented here hold to the ex- 
tent that Sersic's R 1 ' n model provides a good description 
of the various galaxy light-profiles. For early-type galax- 
ies this is known to be the case (e.g., Caon et al. 1993), 
while for late-type galaxies the correction is <0.06 mag 
and one can argue whether or not to apply it. 

For typical elliptical galaxies, Petrosian magnitudes can 
fall shy of the total magnitude by a couple of tenths of a 
magnitude. Application of our method to the SDSS Pet- 
rosian magnitudes yielded total (Sersic) magnitudes accu- 
rate to a couple of hundredths of a magnitude. Moreover, 
it actually provided a more accurate measurement of this 
quantity than obtained from direct Sersic fits to the data 
(Blanton ct al. 2005). 

We have shown how the galaxy concentration can also be 
used to determine the effective half-light radius R c from 
the Petrosian radius Rp, and thereby additionally lead- 
ing to a determination of the associated surface brightness 
terms. We also provide magnitude corrections under the 
assumption that galaxy light-profiles truncate at 5 effec- 
tive radii. 

In regard to the SDSS data, Blanton et al. (2005) showed 



how their Sersic fits to objects with increasingly higher val- 
ues of n resulted in increasingly underestimated effective 
radii. For example, their code obtained radii that were 
~80% of the true effective radii for simulated objects with 
n = 4.5 (i?9o/i?5o = 3.4). Due to this, the concentration- 
corrected effective radii derived from the Petrosian 50% 
radii are in fact preferable to those obtained from the di- 
rect Sersic fits. The same is obviously true for the associ- 
ated surface brightness terms. 

In passing we warn that the use of 3i?p, with l/fy(i?p) = 
0.5, misses more than half of an i? 1 / 4 profile's flux and 
should therefore not be used in the analysis of bright el- 
liptical galaxies. The use of 2i?p, with l/ry(i?p) = 0.2, as 
justified by the SDSS consortium, is suitably appropriate. 
The optimal choice of concentration remains an outstand- 
ing issue. 

While the corrective scheme outlined here is general to 
any choice of Petrosian index and aperture, the tabulated 
values we have provided are specific to two commonly used 
Petrosian systems, namely that used by CAS (Conselice 
2003) and the setup adopted by the SDSS collaboration 
(York et al. 2000). 
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Fig. 1. — Inverted Petrosian index, 1/rj (equation 3), as a function of normalized radius, R/R e , for light— profiles having a range of Scrsic 
shapes n=0.5, 1, 2, 3,... 10 (equation 2). Inset Figure: Number of effective radii that the Petrosian radius Rp — the radius R at which the 
inverted Petrosian index equals some value — corresponds to when the l/r) = 0.2 and 0.5 (i.e., the dot-dashed lines in the main figure). 
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Fig. 2. — Panel a) shows the difference, as a function of light-profile shape 'n', between the Petrosian magnitude (equation 5) inside (i) 
twice the radius where 1/tj=0.2 and (ii) thrice the radius where l/r/=0.5 and the total magnitude obtained by integrating the R 1 '" profile to 
infinity. Panel b) is similar to panel a) but assumes the R}/ n light-profiles are truncated at 5R C . 
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Fig. 3. — Concentration index, defined as the ratio of radii J?go/i?so, is shown as a function of the light-profile shape n. (i) -R90 and R50 
are simply the radii containing 90% and 50% of the total flux (obtained by integrating the Sersic profile to infinity, note that R50 = R c in this 
case), (ii) Similar to (i) except that the total flux is now considered to be only that within 5R e (i.e., the Sersic profile is truncated at 5R C ). 
(iii) -R90 and R50 are now the radii containing 90% and 50% of the flux within NRp , where we have chosen N = 2 and Rp is the Petrosian 
radius such that l/rj{Rp) = 0.2. 
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2 2.5 3 3.5 



R 9o/ R 50 ( <2R p) 

Fig. 4. — a) Magnitude differences shown in Figure 2 for the l/r](Rp) = 0.2 curves versus the associated Petrosian concentration index 
Rgo I R50 within the Petrosian aperture 2Rp (as shown by curve (iii) in Figure 3) . The solid curve gives the correction to obtain the total 
(Sersic) magnitude, while the dash— dot curve gives the correction if the underlying Scrsic profile is truncated at 5R C - b) The ratio of Petrosian 
radii Rp (where l/rj{Rp) = 0.2) to effective radii R c . c) Radius containing 50% of the Petrosian flux, R50, divided by the effective radius 
Rc d) Difference between the surface brightness at R50 and R c . c) Difference between the mean surface brightness within R50 and the mean 
surface brightness within R c . All panels are relative to R90/-R50 within 2Rp and with l/rj(Rp) = 0.2 using the "standard" definition for the 
Petrosian index (equation 3). 
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Fig. 5. — Panel a and b) Difference between the exact magnitude corrections seen in Figure 4a and the approximate magnitude corrections 
given by equation 6. Panel c) Relative error in the radius R x , the approximation used to represent the effective radius Rc (equation 7). Panels 
d-f) Similar to panels a-c) but based upon the SDSS definition of the Petrosian index (equation 10). 



12 



Graham et al. 





0.5 


m 

B 






o 


i 






0.5 




0.5 

















0.5 




2 




1.5 








1 








0.5 




2 




1.5 






\ 


1 


ft;" 






0.5 




1.5 


a 


1 

0.5 









0.5 


a. 


-1 




1.5 


a 


1 

0.5 









A* _ 

a 


0.5 




-1 




1.5 



.. ■ JmU>ami,£uau'..r^t.-... i . - . . ~~~ " *~ " 








a) - 


1 , .1 1— t 1 ( 1 1 1 — 


— 1 1 1 1 1 1 1 1 1 1 






1 . . . 1 


N 


2 2.5 


3 3.5 4 


. 1 1 1 1 1 1 1 


1 , | i i i i . 

c) - 


- ■ - 




1 1 1 1 I I I I 1 1 

/ : 
y : 


' . Ill, 


1 I I I I : 


2 2.5 


3 3.5 4 




1 i | 

e) " 


■■,.]■■. 1 , I I I I I ■ I I 1 X I 




...... . , . ^ ^ 

y 


1 


I 



2.5 3 
^90/^50 ( < 2R P ) 



Fig. 6. — Panel a) The points show the difference between the SDSS Petrosian magnitude mp and the NYU-VAGC Sersic magnitudes mg 
as a function of concentration for 16128 galaxies. The solid curve is the expected difference based on the correction we have formulated to 
obtain the flux beyond the Petrosian aperture. The dashed line reflects the selection boundary imposed on the data by the restriction that 
n < 5.9 in the Sersic fitting process used by Blanton ct al. (2003b). We made a cut at n = 5.8 to avoid any galaxies which may have piled up 
at the upper boundary to n. Panel b) Our corrected Petrosian magnitudes m x versus the Sersic magnitudes. The selection boundary in panel 
a) has been propagated into this panel. Panel c) Petrosian 50— percent radii i?5o divided by the Sersic effective radii R a vs. concentration. 
From Table 1, one can see that the restriction n be less than 5.8 (artificially) prevents R50/R0 from getting smaller than 0.53. Panel d) Our 
corrected values for -R50, denoted here by R x , are shown divided by the Sersic radii R c . Again we have propagated the selection boundary. 
Panels e) and f) are similar to the other panels but show the mean surface brightness. 
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CAS Petrosian: 


Aperture equals 1. 


5i? P (l/r? = 


= 0.2), 77 given in Eq. 3. 
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Note. — Col.(l): Sersic index n. Col. (2): Petrosian radius Rp such that l/ri(Rp) = 0.2. Col. (3): 
Radius containing 50% of the Petrosian flux. Col. (4): Concentration index defined by the ratio of 
radii that contain 90% & 50% (and 80& & 20%) of the Petrosian flux. Col. (5): The inverse of 
column (4). Col. (6): Difference between the surface brightness at R e and the value at i?so- Col. (7): 
Difference between the mean surface brightness inside the radii R c and R50. Col. (8): Difference 
between the Petrosian magnitude and the total magnitude (obtained by integrating the Sersic profile 
to infinity). Col. (9): Difference between the Petrosian magnitude and the Sersic magnitude if the 
nrofile is truncated at 5 7?„_ 



